Finite Drude weight for ID low temperature conductors 
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We apply well established finite temperature Quantum Monte Carlo techniques to one dimensional 
Bose systems with soft and hardcore constraint, as well as to spinless fermion systems. We give 
clear and robust numerical evidence that, as expected, no superfluid density for Bosons or Meissner 
fraction for fermions. is possible at any non zero temperature in one dimensional interacting Bose 
or fermi lattice models, whereas a finite Drude weight is generally observed in gapless systems, in 
partial disagreement to previous expectations. 
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I. INTRODUCTION 

In the last decades there have been a lot of numerical 
and theoretical works to understand the role of strong 
correlation in lattice model Hamiltoniansjii^^^i Re- 
cently this issue has acquired an increasing attention and 
remarkable importance, due to the recent advances in the 
realization of optical lattices. In these experiments ultra- 
cold atoms behave as boson particles trapped on particu- 
lar lattice sites, whereas the interaction and the hopping 
parameters can be tuned continuously. This important 
achievement has opened the possibility to verify directly 
the crucial role played by the electron correlation in very 
important model Hamiltonians defined on a lattice. An 
important example is the realization of a Mott insulating 
state in a system with strong on site repulsion 8,9 . More- 
over quite recently the possibility to include the Fermi 
statistics in optical lattices appears very promising and 
interesting^ 

In ID spinless fermion systems are equivalent to inter- 
acting Bose systems with hard-core constraint and are 
described by the same low energy theory -the Luttinger 
liquid theory-. Indeed this theory holds also for soft- 
core bosons, as shown in Ref.fj]). Therefore, as far as 
the transport properties are concerned one should expect 
the same behavior both for fermions and bosons. On the 
other hand for lattice models, even in absence of disor- 
der, the current does not commute with the Hamiltonian, 
implying its possible decay at finite temperature due to 
the backscattering processes^. In this case the dynam- 
ical current-current correlation function also decays in 
time, leading to a current Fourier transform without S 
function at zero energy, namely without a finite Drude 
weight within the linear response theory. 

Until few decades ago the absence of the Drude weight 
was the expected behavior of all interacting metals in lat- 
tice models or in real solids at finite temperature. How- 
ever a quite clear numerical evidence has been reported 
in Refill that current should not decay in integrable ID 
models, namely for Hamiltonians that can be solved by 
Bethe ansatz techniques in ID. These models essentially 
possess some hidden conservation law, that was conjec- 
tured to forbid the current decay process. 12 ' 13 Later sev- 



eral groups have reproduced this surprising effec t 14 ' 15 , 
with a noticeable exception that a finite Drude weight 
at finite temperature was found also for non-integrable 
models. 15 On the other hand, from purely theoretical 
grounds this issue is not settled yet: in RefQj] it was ar- 
gued that backscattering processes can be effective also 
at finite temperature and in ID non integrable models, 
whereas in RefEE it was proposed that also some par- 
ticular non integrable model could provide a conserved 
current. 

In this work we propose that the general behavior of 
ID gapless systems is eventually characterized by a finite 
Drude weight at finite temperature, and we have found 
no exception in the models that we have studied. This 
conclusion is based on a careful and systematic numerical 
work on fairly generic one dimensional Bose and Fermi 
systems, that all show the same behavior, even though 
strong finite size effects are observed in the non integrable 
cases. 

In the following we investigate the behavior of the 
Drude weight in ID systems in the thermodynamic limit 
and finite temperature. 

Model and Method : We have studied hardcore and 
softcore bosons in a ID lattice with periodic boundary 
conditions. The Hamiltonian studied reads, 



H = f-t(ala i+1 + h.c.) + ^-n^n; 



1) 



+Vriin i+1 + Wriini + 2 - [m, 



(1) 



The sum is over all lattice sites i, aj/ai is the boson 
creation/annihilation operator at site i, henceforth rij is 
the particle number at site i and/i is the chemical poten- 
tial, t is the hopping amplitude which is set to one, U is 
the on-site repulsion, whereas V and W are the nearest 
and the next-nearest neighbor interactions, respectively. 
For hardcore bosons in the U — > oo limit the Hamilto- 
nian can be mapped onto an S = 1/2 spin system with 
Sf = rii — 1/2 and = at. In this work we present 
our results for the half filled case of hardcore and soft- 
core models. Most of our results have been obtained by 
Quantum Monte Carlo (QMC), using the stochastic se- 
ries expansion (SSE)£ii! with the directed loop update^ 8 -. 
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Superfluid density p s (or spin stiffness in the equiva- 
lent spin model), is defined as the second derivative of 
the free energy with respect to a twist in the boundary 
conditions. In order to compute this quantity by QMC, 
it is convenient to apply linear response theory, relating 
this quantity to the current current response function 
A(q,iw„) = / /3 drexp(iw n r)(J(q,r)J(-q,0))/7V, where 
J is the current operator and u> n is Matsubara frequency. 
Then the following expression for the superfluid density 
is obtained: 



p s = (-K) - A(q = 0; ico n = 0) 



(W 2 ) 



(2) 



where (K) is the average kinetic energy per site, oj n = 
2irn//3 arc the Matsubara frequencies and W is the wind- 
ing number. Similarly the Drude weight is obtained with 
the same expression but with a different order in the limit 
to -> and q -> 0, namelyi^^^ 



D = (-K) - RcA( 9 = 0, lo 0). 



(3) 



In SSE one can obtain A very accurately in terms of 
Matsubara frequencies. Therefore analytic continuation 
of the data is required. In order to avoid difficulties of 
extrapolation to iio n — > at large temperatures, we have 
worked at relatively low temperatures ((3 > 10). 

In principle, due to the different order of lim- 
its, the Drude weight and the superfluid density may 
be different when the following quantity remains fi- 
nite in the thermodynamic limit 15 : D — p s = 
Y,E n =E m Pexp(-f3E n )\(ip n \J\ip m )\ 2 /L, where, J is the 
current operator, while E n and \ip n ) are the n th eigen- 
value and eigenstate of the many body system, respec- 
tively. 

The current operator can be written as J(q = 0) = 

iJ2b(H-b ~ Hb) wnere Hb = t a \ a i+i an d b is the bond 
index, corresponding to the site index I. The ensemble 
average of product of two local operators H^ 1 and H^ 2 (r) 
is: 

(H£(t)H£{Q)) = 

1 °° ( T _ ft\n( _ T m \ 

7E E '-(^{H^H^m (A) 



k n,m— 



where r is the imaginary time, Z is partition function 
and the summation over n and m comes from Taylor- 
expansion of e(~^ +T ) H and e~ rH . Following ReftTn the 
relation ([4J can be simplified to 



m=0 



7 { — % - l N b ^^)(5) 
[n s — m — 2y.m\ I w 



where n s is the length of sequence of the local operators 
and it changes in each QMC sampling. N^ 2 ' ai 172 is the 
number of times that two operators H^ 1 and H£ 2 appear 
in this sequence with distance of m local operators, and 



(■■■}w indicates an arithmetic average using configura- 
tions with relative weight W. In this work we introduce 
an efficient way to sample by SSE the current-current 
response function. To this end, we multiply expression 
([5]) by e lu>nT and integrate over the imaginary time r, we 
obtain: 



1 F 1 (m + l,n s ;2i7rn)N^ c 



0\ 



m=0 



(6) 



where 



\_F\{m + 1, n; z) 



D! 



(n — m — 2)!m! J Q 



dx exp(zx)cc m (l — x) 



n — m — 2 



(J) 



is the confluent hypergeometric function. 

Therefore, the current-current correlation acquires 
contributions determined by length of operator string 
n s . All these contributions are stochastically sampled 
in an efficient way, and in each statistical measurement 
the correlation function A(q = 0, ico n ) has the following 
estimator: 



-J a 1 a 2 J2 1 F 1 (m+l,n s ;2i7rn)N^ (8) 



"1 ,cr 2 = ± 



where Nl ia2 = Y\ N blb2 ^ ia2 

Discussion: At zero temperature, for non degenerate 
ground state, the Drude weight and the superfluidity are 
the same. In a ID system at any finite temperature p s is 
expected to be zero in the thermodynamic limit, whereas 
the Drude weight can be non-zero. For hardcore and soft- 
core bosons in a ID lattice, a systematic size scaling of 
the superfluid density p s clearly shows that this quantity 
vanishes in the thermodynamic limit and for any finite 
temperature (see figures [T] and [2]) . Further, we find that, 
for a fixed set of parameters and at half filling, all super- 
fluidity data versus l/L collapse to one curve whenever 
the x-axis is appropriately scaled with the temperature 
T (see figures [1] and [5]) . This analysis suggests the scal- 
ing form p s (P,L) = p s (/3/L). If one takes the order of 
limit T — > after L — > 00, superfluidity remains zero 
even at zero temperature. Notice that by taking first the 
limit T — y and then L — + 00 superfluidity has a finite 
value for the gapless phase, but this is not a signature of 
superfluidity, rather the occurrence of a finite zero tem- 
perature Drude weight. Though in ID is not possible to 
have a finite superfluid density at any non zero tempera- 
ture, several authors have identified the finite zero tem- 
perature Drude weight with the superfluid density for a 
superfluid with vanishing critical temperature. We be- 
lieve that this identification is a bit confusing and there- 
fore we prefer to think about absence of superfluidity and 
superconductivity in ID systems, as commonly reported 
in the textbooks. 

Fig. [3] shows the current-current correlation versus L0 n 
in the metallic and insulating phases of an integrable 
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Hard Core Model, p=\/2, V=1.5 
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FIG. 1: (color online) Superfluid stiffness for an integrable (a) 
and a non- integrable (b) model versus The system size 

L is ranging from 50 to 1200. 



p = 0.5, U = 4t V=0 



■ /S = 75 

* (8=100 
t 0=150 



(8/L 



0.5 



FIG. 2: (color online) Superfluidity of the soft-core bosons 
versus scaled system size at half filling, the on-site interaction 
is U = 4 



model (W — 0, U — oo). The zero-frequency value is 
the superfluid density p s and the limit uj n — > gives the 
Drude weight JJ. For T4 7 = at zero temperature, there 
exists a critical value V c /t — 2 below which the Drude 
weight is finite. In the first case (a) shown in Fig.© 
with V/t = 2 the Drude weight has a finite value at any 
finite temperature, which is consistent with the previous 
worksi^. In the insulating phase (case b) with V/t = 3, 
the superfluid density coincides with the Drude weight 
and they both tend to zero as the system size increases. 

In a non-integrable model such as hard-core bosons 
with nearest and next nearest neighborer interactions 
earlier works have suggested zero Drude weight as system 
size increases. With SSE we can go to very large system 
sizes and low temperatures and check the scaling depen- 
dence of the Drude weight. In Fig. 0] we have plotted 
current-current correlation versus Matsubara frequency 
for different L, and a fixed temperature T = 1/100. As 
shown in the same Figure (JH) we have also found a finite 
Drude weight at finite T in the celebrated Bose-Hubbard 
model with softcore constraint and in several other mod- 
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FIG. 3: (color online) (a) Current-current correlation for an 
integrable model in the metallic phase. The zero frequency 
data shows superfluidity while the extrapolation to n — > is 
the Drude weight. D remains finite with increasing L while 
p s vanishes, (b) In the insulating phase D and p 3 have the 
same value and both tend to zero by increasing L. 
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FIG. 4: (color online) Response function vs. n for (a) hard- 
core bosons with V/t = 1.5, W/t = 1, T/t = 1/100 and (b) 
Bose-Hubbard model with softcore constraint and U/t = 2, 
p/t — —0.4, T/t — 1/25. The system sizes ranges from 
L = 100 to L = 800. 



els (not shown). Although some evidence that few par- 
ticular non integrable models could have a finite Drude 
weight at finite temperature have been reported before, 
here we have found a very convincing evidence that this 
behavior should be generic for ID gapless system regard- 
less from their intcgrability. We have supported this 
statement by state of the art numerical calculations ob- 
tained for very large system sizes and low temperature so 
that all possible extrapolations are perfectly under con- 
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trol. 

In conclusion it turns out that, at low energy, all gap- 
less lattice models studied scale to the Luttinger liquid 
fixed point where the backscattering is a marginally irrel- 
evant coupling and the current is therefore conserved at 
the fixed point. This is therefore a peculiar and generic 
feature of ID. Indeed in 2D systems, such as hardcore 
bosons with n.n. repulsion in a square and triangular 



lattice, we found no difference between p s and D. 
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